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Abstract 

We report the application of the nonlinear a model to study the multi- 
skyrmion problem in the quantum Hall ferromagnet system. We show that 
the ground state of the system can be described by a ferromagnet triangular 
Skyrme lattice near v = 1 where skyrmions are extremely dilute. We find a 
transition into antiferromagnet square lattice by increasing the skyrmion den- 
sity and therefore \v— 1| . We investigate the possibility that the square Skyrme 
lattice deforms to a single skyrmion with the same topological charge when 
the Zeeman energy is extremely smaller than the Coulomb energy. We explic- 
itly show that the energy of a skyrmion with charge two is less than the energy 
of two skyrmions each with charge one when g < g c . By taking the quan- 
tum fluctuations into account, we also argue the possibility of the existence of 
a non-zero temperature Kosterlitz-Thouless and a superconductor-insulator 
phase transition. 
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I. INTRODUCTION 



Recently, there have been a number of experimental f l|-ffl and theoretical in- 
vestigations of the skyrmions in the integer and fractional quantum Hall effect (QHE) in 
2D electron gas (2DEG). It has been shown [H], 11, 24, 27] that the ground state of a two- 
dimensional electron system at the Landau level filling factor near v — 1 is a Skyrme (soliton) 
lattice. Following Sondhi et al. [0] and Moon et al. 0, one may show that the quantum Hall 
ferromagnet can be described by an effective nonlinear a (NLcr) model. In this article we 
also make use of the effective field theory of the quantum Hall NLcr model P,p|,P|-|PT|] to study 
a multi-skyrmion quantum Hall system. Our results for a multi-skyrmion system show that 
the ground state has long range order. We take the advantage of particle-hole symmetry 
II (which is equivalent to skyrmion-antiskyrmion duality in classical field theory) to study 
the skyrmions associated with holes [y < 1). We will show that within V\ c < v ~ 1, the 
ground state is a triangular lattice in order to reduce the Coulomb energy since, the Zeeman 
energy is negligible. Increasing the skyrmion density in such a way that << V2c < v < v\ c , 
leads to the formation of a square lattice. The square lattice is favored by the gradient 
and Zeeman energy in order to allow the antiferromagnetic alignment of the skyrmion's 
orientation via the attractive XY interaction between them. This can be anticipated by 
the frustration of the skyrmionic hedgehog fields within the triangular lattice which is the 
favorable orientation for the Coulomb energy. This leads to a structural phase transition due 
to varying the Landau level filling factor. By taking the quantum fluctuations into account, 
we also argue the possibility of the existence of a non-zero temperature Kosterlitz-Thouless 
and a superconductor-insulator phase transition. We show that the 2DES with the typical 
Zeeman energy at filling factor smaller than one (the validity range of the square lattice) 
are in the insulating phase. The attractive XY interaction makes a tendency for collapsing 
the single skyrmions which are located in different positions. It favors recombination of 
the iV-seperated single skyrmions and forming a single skyrmion with topological charge 
N. The competition between attractive XY interaction and repulsive Coulomb interaction 
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leads to the stability of the multi-skyrmion system in some circumstances. More explicitly, 
we investigate this for a special two single skyrmions system numerically We demonstrate 
that in a special two single-skyrmion case at extremely small values of the Zeeman energy, 
(depends on the filling factor) collapse can occur through the formation of a single skyrmion 
with the topological charge two. 

II. AN OVERVIEW ON NLa-MODEL IN QHF 

In the NLcx-model approach to QHF, the spin of an electron may be described classi- 
cally by a unit vector whose direction may be changed continuously in the space. In this 
representation, the effective Hamiltonian is functional of the unit vector m(r) P|| 

E[m] = Eo[m] + E z [m\ + E c [m] , (1) 

where -Eo[ m ] an d -^[m] are the conventional NLcx model and the Zeeman energy respectively. 
The last term -E c [m] is the Coulomb energy due to the connection between excess electric 
charge density and skyrmion topological density 

£ [m] = ! Jd 2 r(Vm) 2 , (2a) 
E z [m} = ^ F Jd 2 r[l-m z (r)}, (2b) 

EJLM = £jfrJ&4^. (2c) 

Here p s = e 2 /(16v / 27re£o) is the spin stiffness (assuming zero layer thickness for the 2DEG) 
which arises from the Coulomb exchange energy, e is the background dielectric constant of 
the semiconductor, g = ge 2 / (2e£o) is the Zeeman term, g is the effective gyromagnetic ratio, 
and £ is the magnetic length. The charge density is given by p(r) = e a p m(r) ■ 
[d a m(r) x dpm(r)] which is equal to the filling factor times the topological 0(3) spin texture 
density of the quantum Hall ferromagnet (QHF). The total skyrmion charge denoted by Q 
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is an integer- valued topological invariant. It can be determined by the integration upon the 
charge density and classified by homotopy group of a 2D-sphere respectively 

The lowest energy skyrmion solution, m(r), has to satisfy a non-linear differential equa- 
tion which can be obtained by minimizing the energy in Eq.(|I]) with respect to m 

p s (-V 2 + m • V 2 m)m /t - _ ^m,,) - ^-e a p{d a V(r)}(m x dprti)^ = , (3) 

where V(r) is Hartree potential (the exchange potential resides in p s ) 




and p is the skyrmion charge density associated with the minimum energy solution, m(r). 
The solutions of Eq.@ can be classified by the skyrmion charge Q = J dr p(r). From now 
and for the sake of simplicity, we remove the tilde over the classical solution and denote it 
by m. It is easy to find the following equation of motion of the optimal skyrmion by making 
use of cross product of m upon Eq.(|3]) 

d*-J* = 2^g(* x m) A) (5) 

where 

= p s (m x d a m) x - —V(r)e a pdpm x . (6) 

47T 

One may immediately read off from Eq.([5]) that d a = 0. We may define the ground state of 
QHF at v = 1 as a vacuum of skyrmionic spin textures where all spins are aligned along the 
magnetic field direction, i.e. the z-axis. Note that in the absence of the Zeeman energy, the 
alignment of spins along an arbitrary axis occurs due to spontaneous global 0(3) symmetry 
breaking ]6|] hence the minimizing the electrons exchange Coulomb energy. This is also the 
state of spins far from the center of the skyrmions. Therefore the number of skyrmions 
(antiskyrmions) is counted by — N^\. 

Before we consider a lattice of skyrmions we need to have the correct shape of a single 
skyrmion. A single skyrmion is a topological optimal solution of Eq.(^) with Q = 1. For 
our purposes, it is convenient to parameterize the unit vector m by 
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m = (ip x , <p y , yfl-ipip) (7) 

where ip = ip x + i(p y . Near the core of the skyrmion we do not expect that the shape of the 
skyrmion is influenced much by the magnetic field. However, for large distances from the 
core, the Zeeman energy becomes dominant and we need to consider its effect where the 
direction of the unit vector m is close to the vacuum, namely, z-axis. Taking the limit of 
small ip, we can expand the m up to quadratic order in the ip. Putting this in Eq.(0) gives 

m = Jdr (^V 2 ^ + 4^V>) , (8) 
which leads to the equation of motion 
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Introducing k 2 = g / '(2tt£qP s ) we simply have the equation — V 2 ^ + k 2 ^ = 0. We are 
interested in the 'vortex' solution if) = 2d z ^, where d z = (d x + id y )/2 and z = x + iy. (If 
g = 0, this would result in V 2 £ = with the solution £ oc ln(r) and therefore ip oc z/r 2 ). 
Substituting this we find — V 2 £ + k 2 £ = or £ oc e~ Kr / -Jr and therefore that ip ^ ze~ Kr /r 3 ^ 2 
for t — ► oo. The most important of this part is of course the exponential (as opposed to, 
algebraic if g = 0), fall-off PU|. The dynamics of a skyrmion spin texture in NLa-model 



PsV ^ + = 0. (9) 



may be incorporated via the Wess-Zumino action [[R|]. The result of expansion for the single 
skyrmion's Wess-Zumino term is 

Swz = -J^2 J dr J dtm- (d T m x d t m) = -J^ J dt J dr ip(r, t)i^ip(r, t), (10) 

where we keep the quadratic terms in Eq.(|K]). At this level of approximation, the effective 
action may be obtained by Eq.(|) and Eq.(|H]) where the Lagrangian density is 

and S = j dt j drC. The optimal solution of if) is then identical to the solution of the time 
dependent Schrodinger equation where the external potential is proportional to the Zeeman 
splitting factor. Then the single skyrmion behaves approximately like a quantum mechanical 
point particle far from its core. 



III. SKYRMION INTERACTION 



We can now consider the interaction between skyrmions by generalizing our linearized 
energy functional. For that one should note that for the optimal single skyrmion spin texture 
in Eq. (|3|) we may choose a particular orientation. In general the texture can be rotated 
without costing any energy. Since the system shows the U(l) symmetry, any valid skyrmion 
spin texture can be obtained by rotating all spins about the z-axis by angle x- Therefore the 
state, m' = exp(ix/ z )m, is also an optimal solution of the Hamiltonian then I z = —id/dx is 
the generator of the rotation along the z-axis in the internal space. One may expect that it 
contributes to the Hamiltonian through (l/2Ao)(— i<9/<9x — £o) 2 where Ao is the moment of 
inertia of the single skyrmion. This is the leading term of the total energy which is expanded 
with respect to the number of reversed spin, £. In general the dimensional analysis shows 
E z oc £ and E c oc The optimal value of the number of reversed spin, £ = E z /(2g), 

can be evaluated by minimization of the total energy with respect to £ 

6 = ^ / *[1 - ™*( r )] (12) 

where m z (r) is the optimal solution of Eq.(^) |9j corresponding to the given Zeeman splitting 
factor. This leads to the optimal Coulomb and Zeeman energy and then the predicition 
EJE Z = 2 hence A = (d 2 E/dC\ = E z /(6g 2 ). Our goal is now to calculate the interaction 
energy between skyrmions with different orientations. Following Piette et al. ||17|| , we start 
with the conventional NLcr model to find out the proper superposition rule for skyrmions. 
In the absence of the Zeeman and Coulomb energies, the energy functional is scale invariant 
and one may find the optimal solutions analytically. In this case, it can be shown that any 



analytic complex polynomial defines an optimal solution |14]-|16j 

m ^ = (iT^Tp' TT>7p' TTR^) ' (13) 

where w = w x + iw y is any Q-sector analytical function. One may decompose w into a series 
of analytical functions each with Q — 1 
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N 

™ = Y, u r ( 14 ) 
i=i 

Any Q — 1 analytical function represent a single-skyrmion, then the number of skyrmions N 
is clearly the total skyrmions topological charge, i.e. N = Q. Eq.(|T^) denotes a sequence of 
order parameters m.j(r) in configuration space. It is convenient to parameterize the nij(r) 
by 

nij(r) = (smr)j(r) cosQ(r),smr)j(r) smQ(r),cosr]j(r)) , (15) 

where T)j(r) and Cj( r ) are polar and azimuthal field variable associated with jth skyrmion. 
One may define P,|T7| Q = ip — Xj where if is the standard azimuthal angle, e.g. r = 
(r cos <p, r sin <p) and \j is skyrmions orientation and measuring the deviation from the stan- 
dard hedgehog fields. If we have a gas of skyrmions far away from each other, the total 
energy is invariant under variation of skyrmions orientation. For finite separation we expect 
a coupling between skyrmions due to different orientations. Here we consider a situation 
where nij are localized and well separated. This is valid for dilute skyrmions in a quantum 
Hall system, e.g. z/ 2c < v < z/ lc . One may divide the configuration space (R 2 ) into iV regions 
such that Ui is significant in region i and small in others 



1 + Z-VOi 



(16) 



Here \i = (x, y) and 



where j ^ % and $^ = (ip^, <p. 



u jfi = ^ M (18) 



Substituting Eq.([U])-Eq.([yD into Eq . Cp~3|) and expanding up to ipj, we may find m in 
ith region 

m = mj + Qi x nij + -fij x x nxj) + 0(Q 3 ), (19) 



here Eq.(|T9"D describes an infinitesimal rotation of rrij about Qi axis where 

a = im, x {(1 + z ■ naMs - (m ; • $* ff )5}, (20) 

and 

N 

^W = E*iW- (2i) 

Therefore the effect of the other skyrmions on the specific skyrmion is the same as a single 
skyrmion with charge Q — 1 via the effective linear field, $ e fj. In order to study the effect of 
the Zeeman and Coulomb energy, we make the ansatz that the above superposition rule is 
valid for skyrmions far from each other even in the presence of full interaction. This can be 
taken into account by Eq.flTTjp and evaluating the total energy of multi-skyrmion in QHF. 
It is obvious that the total energy can be divided into energies in separated regions. We 
assume that the interaction between the skyrmions is weak, hence fli = m.; x One may 
obtain easily the total energy by redoing the same calculation for all separated regions and 
sum over energies 

E[m] = £ (^N + Jdv ^ ff {-f V 2 + ^h&f) + K ff [m], (22) 

where ip l eS = Q l xeS + i^ l ye s and E^[m] is the effective Coulomb interaction between 
skyrmions. Note that in the absence of the Coulomb term, the saddle point solution as- 
sociated with the scalar field, ip 1 , are vortices, i.e. — V 2 ip l + n 2 ip l = 0. One may divide the 
total energy, Eq . fl22|) , into two parts, the self energy of skyrmions and interaction energy 
which are designated by T and V respectively 

T[m] = £ ffiH + Jdr^i-^V 2 + ^l^'j • (23) 

The first term in Eq.(^) is the total self energy of the isolated skyrmions and the second 
term is the effect of their tail in other regions, i.e., the contribution of their kinetic energy in 
regions far from the core. Here we are interested to study the effect of interaction in a system 
of many skyrmions and their physical relevance. One may find the effective interaction 
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between the skyrmions by making use of the Stokes theorem and neglecting the next nearest 
neighbor terms in Eq. ( |22|) which are described by the terms like j i dr J2k^j ' ^fc and 
h dr Efc^' d a $j ■ d a § k then 

V eS [ m ] = E° cS [ra] + fS[m], (24) 

and 

N 

E° eS H = £ dr d a {J^n iX } = ^ £ [dr^{V 2 - k 2 }#, (25) 
i=i Ji 1 <ij> Ji 

where has been defined for ith skyrmion by Eq.(|6|). i?g ff [m] is the contribution of the 

gradient and Zeeman energy to the effective interaction. It describes a system of interacting 

dipoles. Expanding the charge density of skyrmions, p, in terms of ip, leads to the effective 

Coulomb interaction 



£ e C ff [m] = ^£/^ [dr> 



-,2 



r — r 



+ e ^-p~ E( [drd^(r)d umi ■ Q 3 + f drd^^d^ ■ Q t }, (26) 

where Pi = ^ m« • [da^n-i x dpXJXi]. The first and second term in Eq. (p6|) are the 
Coulomb energy due to the monopole and dipole counterparts of skyrmions respectively. 
The former falls off like R" 1 whereas the latter falls as R~ 2 where the distance between two 
skyrmions is denoted by R. Unlike the dense skyrmions in which the effect of the dipole 
term is crucial, the monopole term dominates for an extremely dilute (i.e. v ~ 1) skyrmions 
system. Neglecting the Coulombic dipole term, one may evaluate the integrals in Eq.(^5|) 
by applying the techniques that were developed for a pair of skyrmions by Piette et al. (17 



Kff[m] = ^Y,fdr [dv' Pf^ + p- 2 £ cos( X , - Xi^MM - R<|), (27) 
Ze 0i^j Ji J i |r-r| 4vr <y> 

where Xj ~ X% an d Rj — Ri describe the in-plane relative orientation and effective distance 

between skyrmion jth and ith. The first term in Eq. ( |2T| ) is the electrostatic monopole term 

independent of the relative orientation. The normalization factor, c = 2y / 2, is calculated 

for a single skyrmion. The second term in Eq. ([27]) describes a classical XY-model where the 
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minimum energy configuration specifies the relative orientation corresponds to Xj ~~ Xi = 11 ■ 
K (x) is the modified Bessel function, hence the coupling between the site i and j decays 
exponentially. We see that there is an exponential decrease of the coupling between the ith 
and jth skyrmions for R ^> The XY term of Eq. fl2?|) favors forming a Q-skyrmion 

(Q = N) by recombining the N single skyrmions, hence the collapse of the lattice. (Note 
that the sign of this term in the antiferromagnet ordering is negative, therefore smaller 
seperation between skyrmions, R, is favorable.) Such a combination costs the Coulomb 
energy. The global minimum solution of the energy functional can be identified by the 
skyrmions self energy and the interaction terms. Apart from the self energy terms which 
leads to an important effect to the recombination of the single skyrmions, one has to take the 
interaction terms into account to predict the proper critical point of this transition. Since 
the Coulomb interaction is long range and the XY term contributes to the total energy as a 
short range interaction, then an antiferromagnet ordering within a square lattice can be the 
global minimum of the total energy, and no combining takes place except for small values of 
the Zeeman energy in agreement with the recent result of Lilliehook et al. ]25|. We examine 
this effect by considering a special case, e.g. two skyrmions system. However, finding the 
global minimume of the total energy for the general iV-skyrmion system as a function of the 
filling factor and the Zeeman splitting factor needs more investigation. An estimate based 
on Eq.(^) shows that the energy of a skyrmion with charge two is lower than the energy 
of two skyrmions each with charge one. Note that the inter skyrmionic XY and Coulomb 
energies are canceled out if g ~ \v — l] 1 ^ 3 at the limit of v — > 1. Let's consider a system 
of two single skyrmions with charge one. The total energy of this system can be obtained 
by employing Eq.(p7[) into Eq.(^). One may compare the total energy of the two single 
skyrmions with the self energy of a skyrmion with charge two in order to estimate the critical 
value of the Zeeman energy 

" 2 _ f ^ f J „/Pi( r )P2(r') c 2 g 



E[m; Q = 2} = 2£[m; Q = 1] + i- / dr [ dr' ^'^y - #o(k|Ri - R*|), (28) 

Zen Ji J2 r — r 47r z 



10 



where E[m; Q] is the self energy of the classical Q-skyrmion. This can be evaluated by 
Eq.(|T]) |2"5| . Results of our numerical calculation for g c for several v are listed in Table |. 
The values of the spin stiffness corresponding to the different fractional filling factors, have 
been evaluated by Moon et al. || numerically using the hypernetic chain approximation. 
We expect that our expression for the inter skyrmionic XY interaction, Eq.(j2"S|), fails for 
filling factors, \v — 1| > Att(k,£o) 2 due to increasing the density of the skyrmions (kR ~ 1). 
At v = 1, the inter skyrmionic interaction is exactly zero and one may find g c = 0.53 x 10 -5 
in agreement with the Lilliehook et al. 25j. At v — 1, the competition between self energies 



yields the transition, since the interaction terms do not play role in this case. However, 
at v < 1 the value of g c is enhanced by the skyrmion interaction. Since the order of the 
Coulomb energy and the short range XY interaction with respect to the number of particles 
are 0(N 2 ) and O(N) respectively, then the effective energy of a mult i- skyrmion system, 
Eq. , becomes dominant by Coulomb energy, hence the stability of the lattice at the 
limit of large N as a local minimum is guaranteed. For larger separations, the XY term of 
Eq.(PTD is therefore negligible compared to the effective Coulomb energy. This implies that 
for a lower density of skyrmions the Coulomb interaction determines the lattice structure and 
leading to a triangular lattice. However, for a more dense system one may expect a different 
stable structural symmetry, in agreement with the calculation of Brey et al. One may 



define a domain of separation for which the square lattice is valid, e.g. z/ 2c < v < v\ c or 
k _1 < R < i? - Here the length scale cut off is denoted by R , beyond that a phase transition 
to a triangular lattice occurs. One may consider the static and dynamical properties of a 



skyrmion lattice [20| using the interaction that we have derived 



V(\Ri -Rj\;xi- Xj) = V Q {\Ri - R,|) + cos(x, - Xj)Vi{\R t - R,|) . (29) 

2 ~ 

In our model V (R) is the electrostatic monopole interaction and V±(R) = ^K (kR). This 



is consistent with the result of variational calculation pO"||M| , p8| . Obviously, the skyrmion 
density may be controlled via the Landau level filling factor. As mentioned above, unlike 
the case where the competition between monopole Coulomb term and the gradient dipole 
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terms is crucial to the stability of the square lattice, the electrostatic monopole interaction is 
the dominant term to determine the structure of the lattice for an extremely dilute [y ~ 1) 
skyrmions. In this situation a specific skyrmion moves in a background which is being 
made by the vacuum of the other skyrmions. Since Ui which has been defined by Eq.([16|), 
is zero every where but in the ith region, then p(r) = J2iPi{ r )- Clearly the total charge 
is the summation upon the individual skyrmion's charge, i.e. Q to t = J2i Q% where Qi is 
the topological charge of localized skyrmion in that region. In this case, the interaction 
between separated skyrmions is independent of the relative orientation. The skyrmions may 
be considered by point particles where Pi(r) = u5(r — Rj) and 

2e , ¥j J J \r - r | Ze i+ . |Kj - ttj\ 

Therefore the ground state is being described by Eq.(^) is clearly a triangular lattice, 
independent of the skyrmion orientation, i.e., a Wigner crystal ||21|| . This is in agreement 
with the result of Green et al. fill. 



One may find the critical filling factor v\ c where a transition into square lattice takes 
place. Note that we have made an approximation associated with R > hT x to demonstrate 
the existence of the square lattice. We may expect the approximation, namely the super- 
position rule of Eq.(^), fails for a high density skyrmions. The crossover between low and 
high density of skyrmions occurs at kR ~ 1 and leading to the lower critical filling factor 
Vic-, i- e - v ic = 1 — g/2ps- Using the typical Zeeman energy, g = 0.015e 2 /eo^o for GaAs and 
taking the advantage of skyrmion-antiskyrmion duality, yields |^2 C — 1| ~ 0.3. Our estimate 
for lower limit of v c shows good agreement with the result of random phase approximation 
of Brey et al. Q. 

The appropriate quantum mechanical spectrum of the Skyrme lattice in our model can 



be obtained through Eq.(27) and the Wess-Zumino action 



N 

E 

i=i 



+ 5>o(|Ri-Ril) 



+ x; ViQRi - Rj\) cos(x, - xj) ) *([R; x];t) = ^n^([R■, *];*), (31) 
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where [R; x] = (R-ij Xi> R-2> X2, Rat; Xat) and = cosifi, RiSmipi). ipi is the standard 
azimuthal angle, indicating the ith skyrmion. A = £o/(3<?) is the moment of inertia of a 
single skyrmion, and £0 = E z /(2g). Obviously, Eq.(|3T]) describes a quantum rotor problem 
for x-degree of freedom ||19|| , consistent with the result of Ref. [[27],[^J. This model has 
been used to describe the superconductor-insulator transition in granular superconductors 
and Josephson junction arrays fL8|. We may expect an insulator-superconductor phase 
transition occurs due to variation of the Zeeman splitting (g) at a given filling factor [y < 1). 
Using the empirical data for practical systems which are consistent with the experiments 
of Bayot et al. |4j (g = 0.015e 2 /e£o ~ 2K at v — 0.8) in which the XY coupling constant, 



J = (2g/n 2 )K (^g/27i£ 2 oPs R) ~ 0.05K, is a fraction of the (2A0)" 1 ~ OAK leads to the 
prediction of the Mott insulating phase for this system where the position of skyrmions 
are fixed. We also expect that the system go through superconducting phase for J greater 
than Aq 1 since the skyrmions earring the electrical charges. The ground state described by 
Eq.(pT|) for skyrmions orientation is a quantum antiferromagnet where the U(l) symmetry 
is broken spontaneously. Here E z = 0.214e 2 /e£o is the Zeeman energy associated with the 
optimal single skyrmion solution corresponding to g = 0.015e 2 /d?o, is evaluated by the 
solving the non-linear, non-local integro-differential equation Eq.(|3|), for a single skyrmion 
numerically. The details of this calculation will be presented elsewhere ||. 

The Hamiltonian, Eq . (|3ll) , may also lead to a finite temperature Kosterlitz-Thouless 
phase transition which might be observable. We may predict a critical temperatures for 
the Kosterlitz-Thouless phase transition based on the classical XY-model. Our calculation 
for the classical 2D nearest neighbor XY-model on the square lattice, shows that the phase 
transition occurs around T c = 0.9J ~ 0.047^ at v = 0.8 (Bq ~ 7T) where we consider the 



square lattice constant, a = 2£ yn/\u — 1|, and the nearest neighbour separation between 
skyrmions, R = a^/2/2. Recent observations of an anomaly in the heat capacity may support 
such a phase transition P|,p6f. However, this interpretation can be controversial since one 
may argue that this effect is a remnant of the nuclear Schottky anomaly JT9] , P? . 
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IV. CONCLUSION 



We showed that for extremely dilute skyrmions, the ground state is a triangular lattice 
to minimize the Coulomb repulsion. However at higher densities, a square lattice forms 
to optimize the spin gradient and Zeeman energies. The ferromagnet triangular lattice 
become frustrated by increasing the skyrmion density at critical filling factor v c and a tran- 
sition occurs into an antiferromagnet square lattice. We also argued the possibility of the 
Kosterlitz-Thouless and superconducting-insulator phase transition. We have shown that 
in some special cases the multi-skyrmionic system becomes unstable with respect to the 
formation of the single-skyrmion system with the same conserved topological charge. 
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TABLES 

TABLE I. The critical 5- factor for several Landau level filling factors are presented. The two 
single-skyrmions system become unstable at g < g c and a transition to a single skyrmion system 
with charge two occurs. 

v p s (e 2 /d ) 9c(e 2 /do) 

1.0 2.49 xl(T 2 0.53 xl0~ 5 

1/3 9.23 xlO" 4 0.19 xl0~ 3 

1/5 2.34 xlO" 4 0.52 xlO" 4 
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